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ABSTRACT 

Context. In the era of large spectroscopic surveys, galactic archaeology aims to understand the formation and evolution of the Milky 
Way by means of large datasets. In particular, the kinematic and chemical study of the thick disc can give valuable information on the 
merging history of the Milky Way. 

Aims. Our aim is to detect and characterise the galactic thick disc chemically and dynamically by analysing F, G and K stars, whose 
atmospheres reflect their initial chemical composition. 

Methods. We performed a spectroscopic survey of nearly 700 stars probing the galactic thick disc far from the solar neighbourhood 
towards the galactic coordinates (/ ~ 277°, b ~ 47°). The derived effective temperatures, surface gravities and overall metallicities 
were then combined with stellar evolution isochrones, radial velocities and proper motions to derive the distances, kinematics and 
orbital parameters of the sample stars. The targets belonging to each galactic component (thin disc, thick disc, halo) were selected 
either on their kinematics or according to their position above the galactic plane, and the vertical gradients were also estimated. 
Results. We present here atmospheric parameters, distances and kinematics for this sample, and a comparison of our kinematic and 
metallicity distributions with the Besanfon model of the Milky Way. The thick disc far from the solar neighbourhood is found to differ 
only slightly from the thick disc properties as derived in the solar vicinity. For regions where the thick disc dominates (1 < Z < 4 kpc), 
we measured vertical velocity and metallicity trends of dV^/dZ = 19 ± 8 km s kpc and d[M/H]/dZ = -0.14 ± 0.05 dex kpc 
respectively. These trends can be explained as a smooth transition between the different galactic components, although intrinsic 
gradients could not be excluded. In addition, a correlation dV^/d[M/H] = -45 ± 12 km s dex between the orbital velocity 
and the metallicity of the thick disc is detected. This gradient is inconsistent with the SDSS photometric survey analysis, which did 
not detect any such trend, and challenges radial migration models of thick disc formation. Estimations of the scale heights and scale 
lengths for different metallicity bins of the thick disc result in consistent values, with hjj ~ 3.4 ± 0.7 kpc, and hz ~ 694 ± 45 pc, 
showing no evidence of relics of destroyed massive satellites. 

Key words. Galaxy: evolution ~ Galaxy: kinematics and dynamics - Stars: abundances ~ Methods: observational 



^ ■ 1. Introduction grate because of resonances with the spiral arms (and the cent ral 

k>( bar), as suggested fo r example bv lSchonrich & BinnevI (l2009l) or 

rS In the paradigm of a dark energy and dark matter dominated iRnSVnr pt al l (l7nfish9 



In the paradigm of a dark energy and dark matter dominated |Ro§kar etall (l2008h ? 
Universe, the process of disc galaxy formation is still very poorly 
understood. For instance, the creation of the galactic thick disc 

still remains a riddle. Was it formed by the accretion of satel- To answer these questions, we would ideally like to tag or 

Utes that h ave deposited thei r debris in a roughly planar con- associate the visible components of the Galaxy to parts of the 

figuration jAbadi et all |2003 >). or was the thin disc heated af- proto-galactic hierarchy All necessary consti-aints can be ob- 

ter successive small accr etions as suggested for example by tained from a detailed analysis of the chemical abundances of 

IVillalobos & Helmi ( 2008)? Did the thick disc stars form in situ, cool and intermediate temperature stars. Indeed, F, G and K type 

through heating of the thin di sc by gas rich mer gers and starburst dwarf stars are particularly useful to study galactic evolution, 

during the merger process (IBrook et al.ll2004 . or did they mi- because they are both numerous and long-Uved, and their at- 
mospheres reflect their initial chemical composition. However, 



* Based on VLT/FLAMES observations collected at the European ^ direct measurement of their spatial distribution requires accu- 

Southern Observatory, proposals 075.B-0610-A & 077.B-0.382. rate estimates of stellar distances, which is a delicate step involv- 

** Tables d El and H are only available in electronic form at the ing (if the parallax is not available) the determination of precise 

CDS via anonymous ftp to cdsarc.u-strasbg.fr (130.79.128.5) or via stellar parameters (effective temperatures, surface gravities and 

|http : //cdsweb . u- strasbg . fr/cgi-bin/qcat?3/A+A/, metal content). 
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That is the reason why most spectroscopic surveys of the 
thick disc up to now are restricted to the solar neighbourhood 
(within about a 500 pc radius) except for those with relatively 
few stars in their samples. Nevertheless, detailed metallicity 
measurements based on such spectroscopic observations of kine- 
matically selected thick disc stars revealed in several studies a 
clear distinction in chemical elements ratios between the metal- 
poor tails of the thi n and the thick discs (Edvardsson et al. 1993; 
Bensbv & Feltzind 120061: iReddv et al.1 12006; Fuhrmann 2008; 
Ruchti et al..,2010t lNavarroetal.ll201lh . Stars belonging to the 
latter are a-enhanced, suggesting a rapid formation of the thick 
disc and a distinct chemical history. Indeed, the [ff/Fe] chem- 
ical index is commonly used to trace the star formation time- 
scale in a system in terms of the distinct roles played by super- 
novae (SNe) of different types in the galactic enrichment. The 
a-elements are produced mainly during Type II SNe explosions 
of massive stars (M > 8M0) on a short time-scale (~ 10^ years), 
whereas iron is also produced by Type la SNe of less massive 
stars on a much longer time-scale (~ lO' years). 

On th e other hand, photometric surveys such as the SDSS 
jYork et al. 2000) explore a much larger volume, but suffer from 
greater uncertainties in the derived parameters. Based on pho- 
tometric metallicities and distanc es of more t h an 2 m illion F/G 
stars up to 8 kpc from the Sun, llvezic et aL I 120081) suggested 
that the transition between the thin and the thick disc can be 
modelled as smooth, vertical shifts of metallicity and velocity 
distributions, challeng ing the view of two dis tinct populations 
that was introduced bv lGilmore & ReidI (Il983l) . 

The goal of this paper is to put additional constraints on the 
vertical properties of the thick disc. For that purpose, we spectro- 
scopically explored the stellar conten ts outside the solar neigh- 
bourhood using the lOiha et al.l d 19961) catalogue. The authors of 
that survey photometrically observed several thousands of stars 
towards the direction of galactic antirotation. They provide the 
Johnson U, B, V magnitudes as well as the proper motions for 
most of the targets up to V~18.5 mag. We selected 689 of them 
based on their magnitudes to probe the galactic thick disc, and 
observed them using the VLT/FLAMES GIRAFFE spectrograph 
in the LR08 setup (covering the infrared ionized calcium triplet). 

The line-of-sight (/ ~ 2 77°, b ~ 47°) was chosen in ac- 
cordance to results found bv lGilmore et al.l (l2002h towards (I ~ 
27 0°, b ~ -45°) and (/ ~ 270°, b ~ 33°) that were confirmed 
bv lWvse et alJ (l2006h towards (/ ~ 260°, b ~ -23°), (/ ~ 104°, 
b ~ 45°), (/ ~ 86°, b ~ 35°), which state that thick disc stars 
farther than 2 kpc from the Sun seem to have a rotational lag 
greater than the canonical thick disc. Indeed, it is commonly ac- 
cepted that the latter lags the Local Standard of Rest (LSR) by 
~ 50 km s"', whereas the authors cited above found a lag twice 
as high (~ 100 km s"') at long distances. Because the angular 
momentum is essentially a conserved quantity in galaxy forma- 
tion, this lagging sub-population has been suggested by the au- 
thors to be the remnants of the last major merger of the Milky 
Way, back to z ~ 2. 

For the analysis of this very substantial sa mple of FLAMES 
spectr a we used the pipeline presented in iKordopatis et alJ 
(1201 11 hereafter paper I). It allowed us to obtain the effective 
temperature (Teff), the surface gravity (log g) and the overall 
metallicit}{3 ([M/H]) for the stars of our sample. They were com- 
bined with the proper motions and the (B-V) colours of the Ojha 
catalogue and our derived radial velocities. Those parameters 



We define the stellar overall metallicity as [M/H]=log(^j^) 



N(H)' 



^°8(jv(Jf) )©' where A' represents the number density and M all elements 
heavier than He . 



were then used to estimate the distances, galactocentric positions 
and kinematics, as well as the orbital parameters (eccentricities, 
apocentric and pericentric distances, angular momenta) of our 
targets. 

The structure of this paper is as follows. The observed stel- 
lar sample is presented in the next section together with the data 
reduction and the radial velocity derivation. The method used 
to determine the atmospheric parameters is reviewed in Sect. [3] 
In Sect. |4]the estimates for the atmospheric parameters, the dis- 
tances and the kinematics are presented. In Sect. |5]a clean sam- 
ple of stars is selected, whose stars are analysed in Sect. |6] with 
respect to their distance above the galactic plane; then w e com- 
pare t hem to the Besanfon model of the Milky Way (RobineTal] 
|2003|) . The galactic components, chosen with a probabilistic ap- 
proach and according to their distance above the galactic plane, 
are also characterised and discussed in relation with thick disc 
formation scenarios. Finally, we estimate in Sect. Q the radial 
scale lengths and scale heights of the thin disc and the thick disc. 

2. Observations, data reduction and radial velocity 
derivation 

The observations were obtained with VLT/FLAMES feeding the 
GIRAFFE spectrograph in MEDUSA mode, that allows a simul- 
taneous allocation of 132 fibres (including sky). The GIRAFFE 
low-resolution grating LR08 (8206-9400 A, R -6500, sam- 
pling=0.2 A) was used during ESO observing periods 75 and 
77 (2005 and 2006, respectively). One of the interesting points 
of that configuration is that it contains the Gaia/RVS wavelength 
range (8475-8745 A), and is similar to its low-resolution mode 
{R -7000). 

In that wavelength range the IR Can triplet (8498.02, 
8542.09, 8662. 14 A) is predominant for most spectral types and 
luminosi ty classes even f or very metal-poor stars (see, for ex- 
ample ,Zwitteretal]|2003)- In addition, these strong features are 
still detectable even at low signal-to-noise ratio (S/N), allowing 
a good radial velocity (Vrad) derivation and an overall metal- 
licity estimation. Paschen lines (for example 8502.5, 8545.4, 
8598.4, 8665.0 A) are visible for stars hotter than G3. The Mg i 
(8807 A) line, which is a useful indicator of surface gravity (see 
Ruck & Smith 1993), is also visible even for low S/N. Finally, 
molecular lines such as TiO and CN can be seen for cooler stars. 

The 689 stars of our sample were selected only from their V- 
magnitudes and the availability of prop er motions measurements 
in the catalogue o f lOihaetalJ (Il996h . They were faint enough 
to probe the galactic thick disc and bright enough to have ac- 
ceptable S/N {my < 18.5 mag). To fill all available fibres of 
MEDUSA, brighter stars in the Ojha catalogue were added to 
our survey (my > 14), resulting in a bimodal distribution in 
magnitudes. We made no colour selection, which resulted in a 
wide range of both log g and T^^fi as we will show in Sect. |3] and 
Sect. ID The magnitude precisions range fro m 0.02 mag for th e 
brightest, to 0.05 mag for the faintest stars dOiha et al.lll996l) . 
Associated errors for the proper motions are estimated to be 
2 mas year"'. 

Eight different fields were observed with two exposures 
each. The spectra were reduced using the Gas gano ESO 
GIRAFFE pipelin^ (version 2.3.0. 'I zzo etal]|2004 with opti- 
mal extraction. The sky-lines were remov ed from our spectra us - 
ing the software developed by M. Irwin ( Battaglia et ani2006l) . 
which is optimal for the wavelength range around the IR Ca n 



|http : //www ■ eso . org/ gasgzino| 
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Fig. 1. Radial velocity uncertainties versus the (B-V) colour for 
the final catalogue (479 stars), as selected in Sect. |5] Dots with 
different colours refer to diflterent metallicity ranges. 



triplet. For each field the four fibres that were allocated to the 
sky were combined to obtain a median sky spectrum and sepa- 
rate the continuum from the sky-line components. Then, the lat- 
ter was cross-correlated to the object spectrum (which was also 
split into continuum and line components) to match the object 
sky-line intensities and positions. Afterwards, the sky was sub- 
tracted by finding the optimal scale factor between the masked 
sky-line and the object spectrum. Individual exposures were then 
cross-correlated using IRAF to put them on the same reference 
frame and to be summed. Before summing them, the cosmic rays 
and the remaining residual sky-lines were removed by compar- 
ing the flux levels of the individual exposures. 

To measure the radial velocities we used a binary mask 
of a KO type star at the LR08 resolution (available from the 
Gene va observatory with the girBLDRS routine, see Rover et al. 
l2002h . The cross-correlation function (CCF) between the ob- 
served spectra and the binary mask was computed, extracting 
from the position of the peak the Vrad value. The spectral fea- 
tures corresponding to the full-width at half-maximum (FWHM) 
of the Can triplet and the Mgi (8807 A) lines were added man- 
ually into this binary mask, because we found that otherwise 
the cross-correlation routine did not converge for the lowest 
S/N spectra (where the spectral lines are hard to identify). The 
wider lines added to the binary mask play a dominant role in the 
cross-correlation method. This implies that no particular con- 
cern has to be raised about errors caused by possible template 
mismatches. The only expected effect is a decrease in the preci- 
sion of the Vrad estimates, owing to a broader CCF. Nevertheless, 
as shown in Paper I, the derivation of the atmospheric parame- 
ters is not altered as long as the uncertainties in Vi^d are less then 
~7-8 km s"' . We found in our survey that the mean error on the 
Vrad is 4.70 km s"' and the standard deviation of the error distri- 
bution is 1 .3 km s"' . In Fig.[T]we plot the estimated errors on the 
radial velocities versus the (B-V) colour of the targets. A colour 
code was added according to the derived metallicity of the stars 



(as found in Sect.[3]l. The largest errors are found for the hotter 
and the most metal-poor stars, as expected, because the mixture 
of broader lines for the former and weaker lines for the latter 
differ the most. The individual values of the heliocentric Vi-ad are 
presented in the online Table [1] and will be used in Sect. l4.3l to 
determine the galactocentric velocities and hence the orbits of 
the stars. 

The spectra were then shifted at the rest frame and linearly 
rebinned by a factor of two, in agreement with Shannon's crite- 
rion. Furthermore, they were resampled to match the sampling 
of the synthetic spectra library we had in our possession, which 
we used to derive the atmospheric parameters (see Paper I). This 
led to a final sampling of 0.4 A and an increased S/N per pixel. 

Then the spectra were cut at the wavelengths 8400-8820 A to 
keep the range with the predominant lines and suffering least by 
CCD spurious effects (border effects, presence of a glow in the 
red part) and possible sky residuals. For this reason, the range be- 
tween 8775 and 8801 A which contains few iron lines compared 
to the possible important sky residuals, was also removed. The 
spectral feature corresponding to the Mg i around ~8807 A was 
nevertheless kept. The cores of the strong Ca ii lines were also 
removed, as explained in Paper I. 

The final spectra, containing 957 pixels, were normalised 
by iteratively fitting a second-order polynomial to the pseudo- 
continuum, asymmetrically clipping off points that lay far from 
the fitting curve. This first normalisation is not very crucial, be- 
cause the pipeline that derives the atmospheric parameters renor- 
malises the spectra with the aid of synthetic spectra (see Paper I). 

Finally, the binary stars or suspected ones, and very low qual- 
ity spectra (owing to a bad location of the fibre on the target) 
were removed at this point from the rest of the sample, lead- 
ing to a final sample of 636 stars. The final measured S/N of 
our spectra varies from ~ 5 to ~ 200 pixeP', with a mean of 
~ 70 pixeP'. Nevertheless, the cumulative smoothing effects of 
spectral interpolation and resampling lead to an over-estimation 
of the S/N by a factor of ~ 1 .4 because of the pixel correlation. 
Hence, the plot of Fig. |2] shows the apparent magnitude versus 
the corrected S/N. 

3. Determination of tlie stellar atmospheric 
parameters 

We used the procedure described in Paper I to derive the atmo- 
spheric parameters of the stars and their associated errors. To 
take into account the overestimation of the S/N owing to the 
pixel resampling, the measured S/N was decreased by a factor 
1.4 at each step of the pipeline. The obtained parameter values 
are shown in the online Table |2l 

We recall that this procedure consists of using two differ- 
ent algorithms simul taneously, MATISS E jRecio-Blanco et all 
2006) and DEGAS (iBiiaoui et alJ l2010h . to iteratively renor- 
malise the spectra and derive the atmospheric parameters of the 
observed targets. The learning phase for the algorithms is based 
on a grid of synthetic spectra covering Teff from 3000 K to 
8000 K, log g from to 5 (cm s^^) and [M/H] from -5 dex 
to -hI.O dex. A coupling between the overall metallicity and the 
a element abundance^ is assumed according to the commonly 
observed enhancements in metal-poor galactic stars: 

- [a/Fel^O.O dex for [M/H] >0.0 dex 

- [a'/Fe]=-0.4x[M/H]dexfor-l <[M/H]<Odex 

^ The chemical species considered as or-elements are O , Ne , Mg , Si , 
S , Ar , Ca and Ti . 
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4.1. Procedure to estimate stellar distances 
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Fig. 2. Apparent magnitude versus measured signal-to-noise 
ratio (per pixel) for all observed sample (689 stars in total). 



- [a/Fe]^+0.4 dex for [M/H] < -1.0 dex. 

At S/N~50 pixer' the pipeline returns typical eiTor^ for 
stars with log g >3.9 and -0.5 <[M/H] < -0.25 dex, of 70 K, 
0.12 dex, 0.09 dex for Teff, log g and [M/H], respectively. For 
more metal -poor stars (-1.5 <[M/H] < -0.5 dex) the achieved 
accuracies are 108 K, 0.17 dex and 0.12 dex. Finally, for stars 
classified as halo giants (T^ff <6000K, log g <3.5, -2.5 <[M/H] 
< -1.25 dex) the typical errors are 94 K, 0.28 dex and 0.17 dex 
(see Paper I). 

This pipeline was successfully tested on the o bserved stel- 
lar libraries of dAUende Prieto et all l2004 and CFLIB 
dValdes et al.l2004l) . showing no particular biases in T^ff or log g. 
Still, a bias of -0.1 dex was found for metal-rich dwarfs, which, 
as discussed in Paper I, we decided not to coiTect because of 
the small range of log g and [M/H] spanned by the libraries. As 
discussed in Sect. 14.11 this possible bias is not expected to intro- 
duce any significant bias in the distance estimates or the derived 
velocities of the present survey. 



4. Determination of the stellar distances, 
kinematics and orbits 

Typical F, G and K main-sequence stars have mv ~ 15- 18 at dis- 
tances of 1-5 kpc, and at least until the ESA/Gaia mission, stel- 
lar distances for targets far from the solar neighbourhood have 
to be determined spectroscopically or photometrically. For in- 
stance, the atmospheric parameters determined in the previous 
section can be projected onto a set of theoretical isochrones to 
derive the absolute magnitudes of the stars. Then we can derive 
the line-of-sight distances using the distance modulus. 



We generated our own set of isochrones by using the YYmix2 
interpolation code, based on the Yonsei-Yale (Y^ ) models (ver- 
sion 2 , [OemarQue et al. 2004) combined with the iLeieune et al.l 
(Il998h colour table. 

We set the youngest age of the isochrones at 2 Gyr, because 
we expect to observe targets from the old thin disc, the thick disc 
and the halo, and very few stars are expected to be younger than 
this. Following Zwitteretal. (2010), we generated isochrones 
with a constant step of 1 Gyr, up to 14 Gyr. In addition, we 
used the full metallicity range of the models, ranging from 
[Fe/H]=-3.0 dex to [Fe/H]=-i-0.8 dex. The step in metallicity is 
constant, equal to 0. 1 dex, smaller than the typical error on the 
derived metallicities of the observed stars. The adopted values 
for the Qf-enhancements at different metallicities are the same as 
those used for the grid of synthetic spectra described in the pre- 
vious section. At the end, a set of 494 isochrones were generated. 

To obtain the ab solute magnitude My, we used the method of 
IZ witter et al.l (12010 ). This procedure consists of finding the most 
likely values of the stellar parameters, given the measured atmo- 
spheric ones, and the time spent by a star in each region of the 
H-R diagram. In practice, we selected the subset of isochrones 
with [M/H]+A[M/H], where A[m/h] is the estimated error on the 
metallicity, for each set of derived Teif, log g and [M/H]. Then 
a Gaussian weight was associated to each point of the selected 
isochrones, which depends on the measured atmospheric param- 
eters and the considered errors (see Eq.lTJ. This criterion allows 
the algorithm to select only the points whose values are close to 
those derived by the pipeline. Still, in low S/N spectra, whose 
parameter errors can be significant, an excessive weight can be 
appointed to rapid evolutionary phases of the isochrones, where 
it is unlikely that there are many stars. To avoid this artefact, we 
associate to the previous weight, another one, that is proportional 
to the tim e spent by a s tar on each part of the H-R diagram. 
Following IZ witter et al.l (1201 Ol) . this weight was set as dm: the 
mass step between two points of the same isochrone. Hence, the 
total weight W for each point on the isochrone can be expressed 
as follows: 



W - dm ■ exp 



-Z 



2A? 



(1) 



defined as the 70% value of the internal error distribution 



where / corresponds to Tetj, log g, [M/H] or (B-V), 0, to the 
parameter values of the points on the isochrones, 0, to the values 
derived from observations (from photometry or spectroscopy) 
and Ag to the associated error of the measurement. The absolute 
magnitude My is then obtained by computing the weighted mean 
of all points of the subset of isochrones. The corresponding 
error is computed by considering the standard deviation of the 
latter. 

To test the adopted procedure, we made exte nsive use of the 
Besangon model of the Milky Way dRobin et al.n2003,) . This is a 
semi-empirical model that includes physical constraints and cur- 
rent knowledge of the formation and evolution of the Galaxy to 
return for a given line-of-sight (los) the atmospheric parameters, 
the positions and the kinematics of simulated st ars in a given 
range of magnitudes. It uses the isochrones of ISchaller et al.l 
(Il992h to simulate the stellar populations of the t hin disc and 
the isochrones of Bergbusc h & VandenbergI (Il992h to simulate 
the thick disc and the halo. Thus, small off'sets may be expected 
when comparing the distances we derived with the Y^ isochrones 
and those deduced from this model. 
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We obtained a simulated catalogue of ~ 4 • lO'' stars towards 
(I ~ 277°, b ~ 47°) in this way and tested the method by com- 
puting the los distances of the targets in two ways: 

1 . The Besan^on atmospheric parameter values were taken and 
associated with errors that the pipeline of Paper I would pro- 
duce for the corresponding stellar type and metallicity at dif- 
ferent S/N. In that way we checked the impact of the weight 
W of Eq. [T] on the distance estimation. We represented in 
the first row of Table |3]the values at 70% of the error distri- 
bution that we found. The results are very encouraging be- 
cause at S/N ~ 20 pixeP' the dispersion of the recovered 
distances with respect to the theoretical ones is only 27%. 
Nevertheless, a bias was found for stars farther than ~8 kpc, 
even at high S/N. For instance, at S/N~ 50 pixeP' the 
distances for the farthest stars are overestimated by ~18%, 
and this offset increases with decreasing S/N (up to 30% at 
S/N~ 10 pixeP'). Indeed, these distant stars are found to be 
metal-poor giants (halo stars), for which the errors in the at- 
mospheric parameters are expected to be larger (e.g.: at S/N~ 
20pixer', the uncertainties are 188 K, 0.57 dex and 0.23 dex 
for Teff, log g and [M/H], respectively). These uncertainties 
will allow a high number of points on the isochrones to have 
roughly the same weight W and to span a relatively wide 
range of magnitudes, because the late-type giant stars have a 
sensitive dependance of absolute magnitude with astrophys- 
ical parameters. This effect will result in a mean My that will 
be lower, designating the stars to a farther distance than they 
should be. This problem is w ell known and has already been 
described in detail in Breddels et al. (2010) and Zwitter et alj 
(2010). As we will see in Sect. 16.2.31 it will concern less than 
4% of our observed sample. 

A disagreement of ~13% towards smaller distances for the 
stars closer than 1 kpc was found. This difference is caused 
by the different sets of isochrones involved, and not by the 
method itself. Indeed, we verified that the absolute magni- 
tudes for main-sequence sta rs at a given metallic ity were al- 
ways higher in the case of ISchaller et alj (Il992h compared 
to the ones. The difference d isappeared for the dwarfs 
ener ated from the isochrones of iBergbusch & Vandenberil 
1992j). 

2. The parameters derived from the spectra of the pseudo- 
stars were taken. For that purpose, we computed degraded 
(S/N~100, 50, 20, 10 pixer') synthetic spectra with the at- 
mospheric parameters of the simulated targets and used the 
automated pipeline to derive their Teff , log g and [M/H] . We 
found that the dispersion on the final distance estimation was 
slightly increased, and that no additional biases were intro- 
duced. Indeed, the projection on the isochrones fixes some 
possible misclassification of the automated spectral classi- 
fication pipeline. For instance, it solves the problem of the 
thickening of the cool part of the main-sequence, as pre- 
sented in Paper I. Therefore, the estimations are quite sim- 
ilar to those obtained from the Besangon values, as seen in 
the last row of Table [3] 

Finally, we tested the effect of a possible metallicity bias on 
the distance determinations. We modified the Besangon metal- 
licities by -0.1 dex to match the suspected bias that the pipeline 
of Paper I exhibited. The errors on the recovered distances were 
only increased by 2% compared to those presented in Table |3] 
We therefore conclude that our distance estimates are not sensi- 
tive to the possible small bias in metallicity that could be present 
in our data. 



Table 3. Values at 70% of the error distribution for the recovered 
distances of the pseudo-stars of the Besangon model. 



S/N (pixer') 


100 


50 


20 


10 


Relative errors (1) 


17% 


20% 


27% 


31% 


Relative errors (2) 


19% 


23% 


35% 


52% 



Notes. Recovered distances (1) using the Besangon atmospheric param- 
eters and (2) using the parameters derived from interpolated synthetic 
spectra at the Besan^on parameters. 
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Fig. 3. Atmospheric stellar parameters for the present sample of 
galactic field stars derived by the pipeline of Paper I and pro- 
jected onto the isochrones. Three isochrones for [M/H]=- 
1.5 dex and 13 Gyr (blue), [M/H]=-0.5 dex and 10 Gyr (green) 
and [M/H]=0 dex and 7 Gyr (red) are represented. The metallic- 
ity values are those derived by the automated pipeline. Typical 
error bars are represented for different regions of the H-R dia- 
gram. 



4.2. Distances of the observed stars 



The observed magnitudes and colours o f lOihaet alJ (119961) were 
dereden ned assuming a me an E(B - V) ~ 0.04 mag and Ay ~ 
0.1 mag (ISchle^el et ay i998). Figures[21|4]and|5]show the atmo- 
spheric parameters projected onto the isochrones, their absolute 
magnitude My, and the histogram of the distances that we ob- 
tained using the method described in the previous section. These 
figures concern the final sample, as selected in Sect.|5] The value 
of Ay should be lower for the closest stars than the adopted one 
because Ay represents the total extinction in the los. However, 
the nearest targets are at Z ~130 pc (D ~ 175 pc), higher than 
most of the dust column, and consequently a single value of Ay is 
a good approximation. In addition, a maximum error of 0. 1 mag 
in Ay will be equivalent to an underestimation of less than 4% in 
the distances, which can be neglected. 
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Fig. 4. Derived absolute magnitude My versus (B-V) colour for 
the present sample of galactic field stars and for the atmospheric 
parameters derived by the pipeline of Paper I and projected on 
the isochrones. The colour codes and the isochrones are the 
same as in Fig. [3] Typical error bars are represented for different 
regions of the H-R diagram. 



Fig. 5. In black we denote the distribution of the line-of-sight 
distances for the subset of the observed stellar sample that meet 
the quality criteria described in Sect. |5] Overplotted in red is 
the distribution of a randomly selected magnitude-limited (my < 
18.5) sample of stars based on the Besan^on catalogue. 



The availability of the distances (D) then allows us 
to estimate their galactic-centred cartesian coordinate^ 

Xcc = Z)cos(Z7)cos(O-X0 (2) 
Ycc = Dcos(b)sm(l) (3) 
Zee = Dsinib), (4) 

where (Xq, Yq,Zq) = (8, 0, 0) kpc (lReidlll993h and /, b are the 
galactic coordinates. 

Table |4] available online, presents the derived values and their 
associated errors (computed analytically, considering no errors 
on / and b). They will be used in Sect. |6]to assign a star to a 
particular galactic component (thin disc, thick disc or halo). 

4.3. Kinematic properties and orbital parameters 

Proper motions combined with distance estimates and radial ve- 
locities provide the information required to calculate the full 
space motions of any star in the Galaxy. For our sample the radial 
velocities were derived from the observed spectra (see Sect. |2]), 
whereas magnitudes, colours and proper motions were taken 
from Ojha et al. (1996). The associated space-velocity compo- 
nents in the galactic cardinal directions, U (towards the galactic 
centre), V (in the direction of the galactic rotation) and W (to- 
wards the north galactic pole) were computed for all the stars 
using the following equations: 



^ We have adopted a right-handed reference frame, with the X axis 
pointed towards the galactic centre. 



U - Vrad cos(b) cos(/) - D/j/, cos(l) sm(b) - Djii sin(/) (5) 
y - Kad cos(b) sin(/) - D^b sin(0 sin(/7) + D/Ji cos(0 (6) 
W = Kad sm(b) + Dfib cos(b), (7) 

where /// and jih are the proper motions on the sk}0. We also 
computed the galactocentric velocities in a cylindrical reference 
frame. In that case, the velocity components are Vr, V^, and Vz, 
with Vz - W. They are computed as follows: 

Vr ^ ^ [Xcc ■ (U + Ue) + Ycc ■ (V + V™, + Vo)] (8) 
= 4 [^GC • iV + Vro, + Vo) - Ycc ■ (U + Uo)] , (9) 



where R = ■\j^cc ^cc planar radial coordinate 

with respect to the galactic centre, (t/f^^, V'®^^, VK^^^) = 
(10.00,5.25,7.17) km s"' is the solar motion decomposed into 
its cardinal directions relative to the Local Standard of Rest 
(LSR, li)e~hneV& BinrrevllT998l) . and Vro, = 220 km s ' is the 
amplitude of the galactic ro tation towards / = 90° and b - 0° 
(lAU 1985 convention; see iKerr & Lvnden-Bell|[T986h . In this 
reference frame a retrograde rotation is indicated by V^ > 
Okms-^ 

The errors on these kinematic data are estimated as follows: 
for each star, we performed 5 10^ Monte-Carlo realisations for 
D, fii, jib and Vrdi, assuming Gaussian distributions around their 
adopted values, with a dispersion according to their estimated 



* One measures /i; = ji'^ cos{b), with yuj the true proper motion in the 
/ direction 
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errors. For every realisation, 6d phase-space parameters and an- 
gular momenta were computed, taking as a final value the mean 
of all realisations, and as an error the standard deviation (see 
Tables |4] and [TJ. The number of realisations was determined to 
achieve a stable result for each measurement. Later on, as de- 
scribed in Sect. 16.21 each of these Monte-Carlo realisations will 
be used to assign for each star a probability of belonging to one 
of the galactic components. 

The median estimated errors on the galactic rotational 
velocities, V^, for metal-rich ([M/H]> -0.4 dex), intermediate 
metallicity (-1 < [M/H] < -0.4 dex) and low-metallicity 
([M/H]< -1 dex) stars are of the order of 10, 18 and 48 km s"', 
respectively. For Vr we reach accuracies of 12, 23 and 
49 km s"', and for Vz we obtain 9, 16 and 41 km s"', respec- 
tively. 



Finally, from the present positions and space-motion vectors 
of the stars one can integrate thek orbits up to several galactic 
revolutions and hence estimate their orbital parameters. Stellar 
orbital eccentricitiefl pericentric and apocentric distances were 
computed. For that purpose we considered a three-component 
Galaxy (dark halo, disc, bulge). The Milky Way potential was 
fixed, and modelled with a Miyamoto-Nagai disc, a Hernq uist 
bulge and a logarithmic dark matter halo (see lHelmill2bo4 and 
references therein for more details). 

Errors on these parameters were propagated from the pre- 
viously found positions and velocities. Typical errors on eccen- 
tricities are ~22%. Results can be seen in Table[T]and in Fig.|6] 
where e has been plotted versus [M/H]. As we will discuss in 
the next section, these eccentricities can be used to identify the 
galac tic components an d distinguish the origins of the thick disc 
stars (ISales et al.1120091) . 

5. Star selection 

We need the cleanest possible sample to confidently describe the 
properties of our observed stars. The final catalogue contains 479 
stars, 452 of which have full 6d phase-space coordinates. Below 
we explain how they have been selected. 

As already mentioned, from the initial 689 targets 53 stars 
that were identified as binaries or had very low quality spectra 
were removed. In addition, we removed the targets for which the 
parameters projected onto the isochrones clearly disagreed with 
the results of the pipeline, according to the following criterion: 



2 2 
^isochrone ^\ 



pipeline 



> 0.25, 



-^ni 



(10) 



pipeline 



where v^. , and , correspond to the differences be- 

^ pipeline ^ isochrone ^ 

tween the observed spectrum and the synthetic templates (re- 
sulting from the pipeline parameters or those projected onto 
the isochrones, respectively). The above condition removed 24 
stars, mainly giants. In addition, 1 1 stars with a projection onto 
an isochrone of 2 Gyr were removed from our survey because 
these young stellar populations are expected neither in the thin 
disc (at those distance above the plane), the thick disc or the 
halo. Removing these targets is justified, because the stars that 
are found to be so young are clearly misclassified by the auto- 
matic pipeline (because of degraded spectral lines, bad normal- 
isation..., or because of the log g-T^ff degeneracy between the 



' We define the eccentricity as e = (rap - rpg)/(rap + rpg) where r^p 
and correspond to the apocentric and pericentric distances of the last 
orbit of each star. 
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Fig. 6. Metallicity versus eccentricity for the 452 stars with full 
6d phase-space coordinates. Typical errors are shown in the up- 
per right corner. 



hot dwarfs and the giants, see Paper I). In addition, this crite- 
rion also removes possible stars on the blue horizontal branch 
(BHB), because the isochrones closest to the BHB will be 
the youngest ones. Finally, we removed 1 1 stars for which the 
spectra had an estimated error in the Vrad > 10 km s"', 43 stars 
with an estimated error on D > 50%, and all remaining stars 
with S/N< 20 pixel '. Indeed, as shown in Paper I, the errors in 
the atmospheric parameters, and hence on the distances are non 
negligible in that case. As a matter of fact, with a selection for 
S/N>20 pixel ', errors less than 190 K, 0.3 dex and 0.2 dex are 
expected for thick disc stars for Teff, log g and [M/H], respec- 
tively, leading to errors smaller than 35% on the distances. 

Figures [3] and |5] show the H-R diagram and the los distance 
distribution for the final catalogue, which contains 479 stars. For 
the 452 stars for which the proper motions were available, Fig.|7] 
illustrates their Vr, V^, and Vz velocity components. The targets 
mainly span los distances from ~175 pc up to ~10 kpc (corre- 
sponding to a distance above the plane Z ~ 7.5 kpc), with only 
1 1 stars reaching up to D ~32 kpc (Z ~ 24 kpc). As described 
in Sect. 14.11 these very distant stars are thought to have strongly 
overestimated distances (up to 30%), and hence one should be 
careful concerning the conclusions obtained for them. In addi- 
tion, the excess of stars seen in the Vz versus Vr plot of Fig. |7] 
towards negative Vz and Vr is caused by small number statistics 
rather than a stellar stream. 

The metallicity clearly decreases with increasing Z (Fig.[8]l, 
as expected, because the proportion of thin disc, thick disc and 
halo stars changes. Nevertheless, we can notice in Fig.|8]a lack 
of metal-poor stars ([M/H]< -1.5 dex) at heights greater than 
~3 kpc from the plane. One could argue that these missing 
stars are either misclassified, unobserved, or removed because 
of one of the quality criteria cited previously. The analysis of 
the magnitudes and the Vrad of the removed stars (empty cir- 
cles in Fig. |9]l shows that these targets are mainly faint stars 
(m^ > 18), with typical velocity values corresponding to the 
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Fig. 7. Galactocentric velocities computed for the observed sam- 
ple as defined in Sect. |5] Typical errors are represented in the 
upper right corners. The colour code according to the metallicity 
is the same as in Fig|3] 



Fig. 8. Metallicity versus distance above the plane (top), versus 
heliocentric radial velocity (y,ad, middle plot) and versus rota- 
tional velocity (V^, lower plot), for the considered sample. Mean 
estimated errors are represented for three different heights or 
metallicity ranges. 
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Fig. 9. Apparent magnitude versus derived heliocentric radial ve- 
locity (y,ad) of the entire sample. Empty circles represent the re- 
jected targets, whereas the filled ones represent the selected stars 
for the final catalogue (see Sect. |5]for more details). Most of the 
faint targets with Vi^d ^ 180 km s"' (typically, the halo stars) 
were removed. 

halo (Vrad > 180 km s"'). This is strong evidence suggesting that 
some metal-poor stars have been observed, but were removed af- 
terwards from the final catalogue. Indeed, the poor accuracies 
achieved for this type of stars (spectral signatures lost in the 
noise, see Paper I) lead to their removal when selecting a clean 
sample for the present analysis. 

Finally, the y^d distribution of the targets extends up to 
Vjjjj = -1-400 km s"' (see Fig. [8]). At these galactic latitudes 
(b ~ 47°) Vrad ^ 300 km s"' corresponds to halo stars with retro- 
grade orbits. The low metallicities of these stars combined with 
the derived V^-velocity component confirm this. These stars will 
be discussed in more detail in Sect. 16.2.31 

6. Characterisation of the observed stellar sample 

To compare our observations with galactic models, we created a 
catalogue of pseudo-stars with the same properties as ours. For 
that purpose we obtained a complete (up to my-\9, mB-2Q) cat- 
alogue of 8 ■ lO'* simulated stars towards our los from the website 
of the Besangon galactic mode0 dRobin et al.l i2003). The latter 
supposes a thick disc formed through one or successive merger 
processes, resulting in a unique age population (11 Gyr) with a 
scale height of 800 pc, a local density of 6.2%, and no vertical 
gradients in metallicity or rotational velocity. 

To run the model, we considered a mean diffuse absorp- 
tion of 0.7 mag kpc We applied parabolic photometric er- 
rors as a function of the magnitude, and a mean error on the 
proper motions of 2 mas year ', according to the values given 
by Ojhaet al. (1996). The input errors for the Vrad were those 
derived in Sect. |2] We selected stars with a Monte-Carlo routine 
to obtain the same my and (B-V) distributions as in the FLAMES 



survey and the same ratio between giant and main-sequence stars 
(resulting from our final catalogue). This strict selection finally 
kept ~ 4 ■ 10^ simulated Besangon stars (called raw Besangon 
catalogue hereafter). Below we will first use this mock catalogue 
to help us in interpreting the vertical properties of the observed 
data (Sect. l6Tt . Then, according to the hints given by the verti- 
cal study, the stars belonging to each galactic component will be 
selected (Sect. 16.2b . which will lead to a characterisation of the 
thin disc, thick disc and halo. 

We recall that the model returns the cartesian velocities U, V, 
W, and to obtain the values Vr, in the cylindrical frame, one 
must use Eqs. [8]and|9l Both of them suppose a solar velocity 
and a galactocentric distance for the Sun. The values adopted by 
the Besan§on model are not the commonly admitted ones (Rq = 
8.5 kpc, Ue = 10.3kms-i, Vg = 6.3 km s"'. We = 5.9 km s"'), 
and for that reason, we preferred to compare only the cartesian 
velocities and not the best-suited cylindrical frame. 

6. 1 . Study according to the distance above tlie galactic plane 

First, we divided the observed sample into stars lying closer and 
farther than 1 kpc from the galactic plane, with 201 and 251 tar- 
gets, respectively. We recall that this cut at Z = 1 k pc roughly 
corresponds to the scale height of the thick disc dJuric et al 



,http : //model . obs-besancon. fr| 



2008) and to the threshold height above which Gilmore et aL 
(2002) identified a lagging population towards lower galactic lat- 
itudes (b ~ 33°). 

The left-side histograms of Fig. [TOl represent a comparison 
between the model and the observations for metallicity, rota- 
tional velocity and Viad of the closest stars. They show that the 
model reproduces the observations close to the plane fairly well, 
except for the position of the metallicity peak, for which we find 
a more metal -poor distribution, by 0. 15 dex. 

On the other hand, for the more distant set (right-side his- 
tograms), several adjustments are needed to bring the metal- 
licities and the velocities of the model to a better agreement 
with the observations. Indeed, the model considers a thick disc 
and a halo following Gaussian metallicity distributions, with a 
mean [mIh]td = -0.78 dex, and [M/H]^ = -1.78 dex, re- 
spectively. These values are too metal-poor compared to the 
mean thic k disc and mean halo metallicities found in the liter- 
ature (see ISoubiran etaTI 120031: iFuhrmannlllool ICarollo et all 
I2010h . We find that a thick disc and a halo more metal-rich than 
the model default values by 0.3 dex and 0.2 dex respectively, 
lead to a much better agreement of the metallicity distribution, 
as can be seen from the comparison between the red dashed 
histogram (the observations) and the black one (the modified 
model) in the first row plots of Fig. [10] We therefore propose 
[M/H]td = -0.48 dex, and [M/H]„ = -1.58 dex. We stress 
that the proposed value for the halo can suffer from the removal 
of some of the metal-poor stars, as described in Sect. |5] 

As far as the velocity distributions are concerned, the model 
still nicely represents the U and W distributions far from the 
plane, but not the velocity in the direction of galactic rotation 
(V). Indeed, we find a shifted distribution with a peak around 
-70 km s ', lower by 20 km s ' compared to the predicted one. 
This ~ 20 km s"' difference is also seen for V^, suggesting that 
its origin is not a spurious effect caused by the local, cartesian 
reference frame. Nevertheless, this effect is hardly visible from 
the radial velocity distribution, and the lag is much less clear 
compared to that visible in Fig. 2 of Gilm ore et al. (2002). This 
might be partly because at higher latitudes the contribution of V 
to the Vi-ad is less important. 
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Fig. 10. Comparisons between our observations (red dashed histograms) and the Besan^on model predictions (black solid his- 
tograms) for the metallicities, rotational (V in cartesian and in cylindrical coordinates system) and radial velocities. The left-side 
plots include the stars that are lower than 1 kpc, whereas the right-side plots include the stars farther than 1 kpc from the plane. 
Dotted black histogram corresponds to the raw Besan§on model (as downloaded from the web, but biased to match our magni- 
tude distribution), whereas the black continuous line corresponds to a model with a richer by 0.3 and 0.2 dex thick disc and halo, 
respectively. 
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To better characterise the observed sample, we studied the 
metallicity and the velocity distributions for naiTower height bins 
up to 4 kpc from the plane. The size of the bins was chosen to 
include at least 20 stars. To fully take into account the uncertain- 
ties in the computed positions and metallicities of the observed 
stars, the plots of Fig. [TT] were obtained from 5 ■ 10^^ Monte- 
Carlo realisations on both parameters {i.e. D and [M/H]). For 
each realisation we computed the new velocities and measured 
the median metallicity and the median V of the stars inside each 
bin. The mean value and the Icr uncertainties rising from the 
Monte-Carlo realisations were plotted. Finally, cr,, was obtained 
by computing the mean robust deviation of V inside each bin. 

A clear change of regime is observed in Fig.[TT]for [M/H], 
V and CTy around 1 kpc, i.e. where the transition between the 
thin and the thick disc is expected. Closer than 1 kpc the metal- 
licity trend is flat, whereas for the more distant stars a gradient 
d[M/H]/dZ = -0.14 ± 0.05 dex kpc is measured This gra- 
dient, agre es (within the erro rs) with that found by iKatz et all 
(1201 lb and lRuchti et all (1201 Ih . Yet this gradient might be under- 
estimated because as noticed previously, our sample might be 
lacking low-metallicity stars at long distances. Velocity gradi- 
ents of dV^/dZ = 19 ± 8 km s kpc and davJdZ = 
9 ± 7 km s~'kpc ' are observed for the stars between 1 and 
4 kpc, which agrees r elativ ely well with the result s found by 
ICasetti-Dinescu et all (1201 lb and lGirard etal] (|2006|) . 

To better understand the origin of these trends, we used 
once more the Besan9on model and compared it to our obser- 
vations. As suspected from the study above and below 1 kpc, 
the raw Besan^on model does not mimic our data correctly 
(Fig. [TT] red diamonds). On the other hand, a more metal-rich 
thick disc and halo, and a thick disc that lags behind the LSR 
by y = -70 km s"', as suggested previously, leads to a much 
better agreement between the observations and the predictions 
(blue triangles). In particular, the measured trends seem to be 
explained in our mock model as a smooth transition between the 
different galactic populations. 

Nevertheless, even with these adjustments the mock model 
does not represent all results we obtained. In particular, bins 
around Z = 1 kpc still clearly disagree, corresponding to the 
height where the thick disc is expected to become the domi- 
nant population. In addition, the plateau at the high-metallicity 
regime is not well modelled either, suggesting a local density of 
the thick disc higher than the one assumed for the model, cor- 
related perhaps with a different scale height. Indeed, the model 
essentially predicts the number of thick disc stars seen above 
one scale height, and thus the adopted thick disc scale height 
and Z = normalisation are degenerate. 

Moreover, a clear correlation between and [M/H] is found 
for the stars between 0.8 and 2 kpc {dV ^/d[M/H] ^ -45 + 
12 km s ~' dex ~') , in agreement with ISpagna et al.l (1201 Ol) 
and iLee etal. ( 20111) ) , thoti gh in disagreement with the SDSS 
view of [ivezic et al.l ('2008'), based on photometric metallici- 
ties. We recall that according to the radial migration scenarios 
dRoskar et al.1 l2008t (Schonrich & BinnevI 12009) . no or only a 
very small correlation is expected in the transition region, be- 
cause the older stars that compose the thick disc have been radi- 
ally well mixed. 

6.2. Characterisation of tiie galactic components 

We now aim to select the thin disc, thick disc and halo members 
to characterise these three galactic components. We recall that 
there is no obvious predetermined way to define a sample of 
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Fig. 11. Median metallicities and V-velocities at diff'erent height 
bins above the galactic plane. Each point has Icr uncertainties, 
obtained from 5 10^ Monte-Carlo realisations on the position 
of the stars (and hence their velocities) and their metallicities. 
Red diamonds represent the predictions of the raw Besan9on 
model, whereas the blue triangles represent a modified model 
with metal-richer thick disc and halo, and a thick disc that lags 
behind the LSR by V = -70 km s ^. 
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Table 5. Adopted values for the characterisation of the galactic 
components 



Galactic 




CTv 






component 


(kms-') 


(km s"') 


(kms-') 


(kms-') 


Thin disc 


35 


20 


16 


-15 


Thick disc 


67 


38 


35 


-70 


Halo 


160 


90 


90 


-220 



Table 6. Galactic population identification obtained from syn- 
thetic spectra of the Besan§on sample simulated with a signal- 
to-noise ratio of approximately 20 pixel"' 



thin disc 


thick disc 


halo 


(1) 


(2) 


(1) 


(2) 


(1) 


(2) 


84% D 


77% D 


25% D 


15% D 


1% D 


0%D 


16% TD 


23% TD 


73% TD 


78% TD 


56% TD 


34% TD 


0% H 


0%H 


2% H 


7%H 


43% H 


66% H 



Notes. (1) Selection according to velocities of the stars (probabilistic 
approach). (2) Selection according to Z-distance. D, TD and H, corre- 
spond to thin disc, thick disc and halo stars, respectively, as identified 
by the Besanfon model. 



a purely single galactic component. Any attempt will produce 
samples contaminated by the other populations. 

W e used the kin ematic app roach of Soubiran & GirardI 
(12005): Bens bv & Feltz ing ( 20061): iRuchti et al.l(l2010h to select 
the stars belonging to each galactic component (hereafter also 
called probabilistic approach). We recall that the advantage of 
this method is that no scale height is assumed for any popula- 
tion. The inconvenience is that it assumes a velocity ellipsoid, 
forcing in that way to find member candidates within the given 
dispersions (introducing for example biases in the eccentricity 
distributions because the cold stars like thin disc's ones will 
have low eccentricities and very hot stars vey high eccentrici- 
ties). This may lead to systematic misclassifications where the 
velocity distributions overlap, or if the assumed distributions are 
not the expected ones {i.e. non Gaussian or with significantly 
different means and/or dispersions). 

In practice, for each set of U, V and W a membership prob- 
ability is computed according to the following equation: 



P = 



1 



exp 



(V - Vlag) 

2ctI 



w_ 

2^1 



2 \ 



(11) 



The adopted values of ctu, cry, ctw and Viag are the ones rep- 
resented in Table |5] They are taken from Bensby & Feltzingj 
(I2OO6I) . except for the Viag of the thick disc, in which case we 
used the value suggested in Sect. 16. II of -70 km s The prob- 
ability of belonging to one of the components has to be signifi- 
cantly higher than the probability of belonging to the others, to 
assign a target to it. In addition, one has to take into account 
the uncertainties on the measured velocities, to avoid massive 
misclassifications. In practice, instead of assigning a component 
based on the final values of U, V and W (presented in the online 
Table HJ, we performed it for each of the 5 • 10^ Monte-Carlo 
realisations used to compute these values (see Sect. 14.3b . The 
adopted threshold in probability ratio above which a realisation 
is assigned to a component was fixed to four (we checked that 
our results are not alfected by the adopted threshold, see below). 
The final assignment was then made by selecting the component 
which occurred in most realisations. 



We tested this approach with the parameters obtained from 
the synthetic spectra of the Besangon sample, presented in 
Sect. 14.11 The adopted velocity ellipsoids are the same as in 
Table |5] but in that particular case, we chose a Viag compat- 
ible with the model, of -46 kms Results for spectra with 
S/N~20 pixel ' are shown in Table|6] where we checked the true 
membership of each assigned star (labelled (1)). This approach 
identifies the thin disc targets quite well (first column), but 25% 
of the thick disc candidates are actually thin disc members (third 
column). Furthermore, the majority of the identified halo stars 
are, in fact, thick disc members. 

The assignment to a component for most of the stars is inde- 
pendent of the adopted probability threshold (in our case: four). 
Indeed, for the majority of the stars more than half of the Monte- 
Carlo realisations result in the same assignment. Increasing 
the probability threshold will only increase the realisations for 
which the assignment will not be obtained. Hence, the compo- 
nent which will have occurred in the majority of realisations will 
remain unchanged, keeping that way the same membership as- 
signment. Finally, we stress that the results ob tained when re- 
placing the velocity ellipsoids of iBensbv & Feltzing. (200 6) with 
those of ISoubiran et al.1 ( l2003l) or lCarollo et al.1 (l2010h" remained 
identical. Indeed, the differences in the velocity ellipsoids or the 
rotational lag are not dilferent enough (few km s ') to introduce 
changes in the final candidate selection. 

To obtain a sample of thick disc and halo stars as pure as pos- 
sible, we found that a candidate selection based on the distance 
above the galactic plane was preferred (like in Di erickx et al.l 
2010, for example). In this case, one has to adopt some val- 
ues for the scale heights of each population, and a normalisa- 
tion factor to estimate the pollution of the other components. 
Scale heights are a matter of debate, especially for the thick disc 
and the halo, but it is generally admitted that stars lying farther 
than -1-2 kpc and closer t han ~4-5 kpc from the galactic plane 
are thick disc dominated (ISiegel etal.ll2002t lJurid et aLlfeoOSl: 
Ide Jong etal.1 120101) . Results obtained for the same Besangon 
sample with this distance selection are also shown in Table |6] 
under the label (2). In that case, only 15% of the thick disc stars 
are misclassified as thin disc members. In addition, the ratio of 
recovered halo targets has increased from 43% to 66% compared 
to the probabilistic approach. 

Keeping these results in mind, we decided below, to investi- 
gate the results obtained by both methods (Z selection and proba- 
bilistic approach) for our observed data. We performed Gaussian 
fits of the distributions of U, V, W and [M/H] for each galactic 
component, and we discuss the mean values and dispersions for 
each method. Results and their Icr uncertainties are shown in 
Table E] 

The probabilistic approach assigned 154 stars to the thin 
disc, 193 stars to the thick disc and 105 stars to the halo. On 
the other hand, in the Z selection the 163 stars lying closer than 
Z = 800 pc were assigned to the thin disc, the 187 stars lying 
between 1 < Z < 4 kpc to the thick disc, and the 45 stars above 
5 kpc from the plane to the halo. 

6.2.1. The thin disc 

The values found with the kinematic approach shown in 
Tabl e [7] are consistent with the properties o f the old thin 
disc (IVallenari et al.ll2006t ISoubiran et al.ll2003l and references 
therein). The mean eccentricity for the thin disc stars is 0.14, 
with a dispersion of 0.06 (Fig.fTJb. though we remind the reader 
that these values are expected to be greater in reality. Indeed, the 
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high eccentricity tail of the thin disc cannot be selected with the 
adopted procedure because of the biases introduced to the orbital 
parameters by the kinematic selection. 

The most distant thin disc star is found at Z~ 1874 ± 224 pc. 
We found one thin disc candidate with [M/H]~ -1.5 dex, at 
Z ~ 250 pc. This star has Vnd - 1-3 ± 4.5 kms"', which 
is typical of the thin disc. Nevertheless, we rather suspect 
this target to be thick disc member, whose kinematics occupy 
the wings of the thick disc distribution function in a region 
that overlaps the thin disc distribution function. This nicely 
illustrates the limitations of our adopted Gaussian model for 
the distribution functions. A measurement of its a-elements 
abundances would increase the dimensionality available for 
population classification, and so might help to identify its 
membership with confidence. 

A selection according to Z-distance results in a hotter veloc- 
ity ellipsoid and a lower mean metallicity compared to the kine- 
matic selection. The higher lag and the lower metallicity found 
with this method is probably caused by the pollution from the 
thick disc. This is also suggested from the eccentricity distri- 
bution of Fig. [141 where one can notice that the thin disc has 
an anomalously high number of stars with e > 0.2. Indeed, we 
recall that in that case, stars up to Z = 800 pc are considered 
to be thin disc members. This model is of course dynamically 
unphysical, and is adopted here merely as a limiting and con- 
venient case to illustrate different population classification out- 
comes. The contamination from the thick disc will thus depend 
mainly on the local density of the latter, whose values are found 
to va ry in the literature from 2% up to 12% (see lArnadottir et al] 
12009 , for a review of the normalisation factors). 

Considering a thin disc with = -211 kms"' (as found 
with the kinematic approach) and a thick disc with - 
-166 km s ', the contamination caused by the latter should be 
~19% to recover the lag of ~ -15 kms"' that we measure. 
Roughly the same result (18%) is obtained when looking for the 
amount of thick disc contaminators which is needed to pass from 
[M/H] = -0.22 dex (kinematic approach) to -0.27 dex (Z se- 
lection) with a canonical thick disc metallicity of -0.5 dex. 

Because the kinematics are reasonably well established in 
the solar neighbourhood, we conclude for the thin disc that the 
probabilistic approach should return more robust results than the 
Z-distance selection. 



6.2.2. The thick disc 

The mean V-velocity (V = -63 + 2 kms"', - -166 + 
2 km s see Tables|7]and[8]l found with the kinematic approach 
is slightly higher than the typical canonical thic k disc value, 
which less than -50 km s ' (see, for example, IWilson et alJ 
12010 ). The mean eccentricity is 0.33, with a dispersion of 0.13 
(see Fig. [T3] ). The closest thick disc star is found at Z = 
129 ± 10 pc, and the most distant one at Z~ 5.33 + 1.17 kpc. 
We find that the metallicity of the thick disc extends from ~ 
-1.8 + 0.1 dex up to super-solar values of ~ +0.25 ±0.1 dex, 
and that [M/H] - -0.41 + 0.02 dex. Let us note again that the 
mean metallicity might be overestimated owing to the contribu- 
tion of the thin disc (misclassifications for the stars lying in the 
low velocity tails), and that the kinematic selection introduces 
biases in the eccentricity distribution. 

On the other hand, the Z selection results in Vr, V^, Vz veloc- 
ity distributions which are compatible with those found proba- 
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Fig. 12. V-velocity dispersion for different metallicity bins for 
the stars lying between 1 < Z < 4 kpc above the plane. Red 
diamonds represent the raw Besan9on model, whereas the blue 
triangles represent the model with our preferred model, which 
has a mean metallicity -0.5 dex, and a mean galactic rotational 
velocity V = -70 km s"' for the thick disc. 

bilisticallj0. Nevertheless, a slightly lower metallicity ([M/H] - 
-0.45 dex) is found, which, interestingly, agrees well with the 
results suggested in Sect. 16.11 where we compared our results 
with those of the Besancon model. 

The argument that the thick disc is a distinct population 
compared to the thin disc is even more accredited by the plot 
of (Tv versus [M/H] (Fig. [T2l i. obtained for the stars lying be- 
tween 1 and 4 kpc. We see that for typical thick disc metallici- 
ties (-1 <[M/H]< -0.2 dex) the velocity dispersion is constant 
within the errors, around cry ~55 km s"'. At these heights and 
metallicities the pollution from the other components is expected 
to be very low (too high and too metal-poor for the thin disc, too 
low and too metal-rich for the halo). Hence, this plot seems to 
highlight the intrinsic velocity dispersion of the thick disc. But 
we cannot rule out the existence of an intrinsic vertical gradient 
in the metallicity and rotational velocity of the thick disc. Higher 
resolution spectroscopy and/or higher number of targets like the 
forthcoming Gaia-ESO survey will help to answer this question, 
though. A cleaner thick disc sample might be obtained by sepa- 
rating the thin disc from the thick disc based on the [a/Fe] ratio, 
and the halo from the thick disc by means of higher statistics. 

Finally, the orbital eccentricities of the thick disc obtained 
with the two methods have to be discussed (see Fig. [T3] and 
Fig. [T4b . Very many authors have recognised that the shape of 
the eccentricity distribution can give valuable information about 
the formation mechanism of the thick disc an d the accreted or 
not accreted origin of its stars (Sales et al. 2009: lDrMatteo et al.l 
1201 ll) . Although for both selection methods the distribution tails 



' We note though that we found non-zero values for and Vz- We 
searched for possible reasons for these offsets without success. These 
mean values could be possibly caused by small zero-point errors on the 
proper motions. 
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can suffer from individual uncertainties or wrong population as- 
signments, the position of the peak is reliable information ac- 
cording to these authors. The eccentricities in Fig.[T3]suffer from 
selection biases, as discussed in the previous sections, because 
the membership was decided on kinematic criteria. Hence, no 
conclusions should be taken based on these results. 

On the other hand, Fig.[T4]suggests a peak at rather low val- 
ues, around e ~ .3. Hence, the pure accretion scenario, as in 
lAbadi et all ( l2003h . which supposes a peak around e ~ 0.5, can 
be ruled out (we note thoug h that these results have been recently 
debated by iNavarr 0I I20T 1|). Nevertheless, no separation can be 
made between the heating, migration or merger scenarios, which 
differ mainly in the detailed shape, especially in the breadth of 
the distribution. 

We note, however, that Fig. [14] suggests a broader distribu- 
tion than has been found so far for more local samples. This is 
very interesting because it might suggest that accretion has con- 
tributed slightly more compared to local sample, and that the 
dominant mechanism for the formation of the thick disk may 
well have varied with distance from the galactic center 

6.2.3. The halo 

The addition of the individual errors on Vrad, pth and D leads 
to an unreliable selection of halo stars based only on their kine- 
matics. Indeed, we recall that our tests conducted on synthetic 
spectra of the mock Besangon catalogue showed that a proba- 
bilistic assignment of halo membership inevitably generates a 
sample strongly polluted by the thick disc. 

To minimise this pollution, our tests suggest to perform a se- 
lection on Z instead. For the 45 stars being more distant from 
the plane than Z - 5 kpc, the mean eccentricity is 0.69, with 
a dispersion of 0.16. We find a mean metallicity of {MjH} - 
-0.92 + 0.06 dex, with ctimih] - 0.76 + 0.06 dex. These values 
are too metal-rich, compared to t he commonly admi tted value 
of ~ -1.6 dex for the inner halo dCarollo et al.ll20T^ . In addi- 
tion, we find mean Vz velocities significantly different from zero, 
which is not expected. These odd values are likely caused by the 
large uncertainty expected for the most distant stars for both the 
distance and hence the velocities to our star selection criterion, 
which created a biased and unreliable halo sample (see Fig|9j 
and probably to an underestimation of the errors for the proper 
motions for the most distant stars (we recall that we assumed a 
constant eiTor of 2 mas/year). 

This could also explain the fairly large velocity dispersions 
that we find compared to those expected for the halo. For ex- 
ample, using SEGUE spectra, CaroUo et al. (2010) found for the 
inner halo (cry, cry, (Tw) = (150 + 2, 95 + 2, 85 + 1). 

We note, however, that another possible solution to explain 
a part of this discrepancy can be found in Mizutani et al. (2003 
who suggested that the retrograde stars in iGilmore et al.l (1200' 
could be the debris of u Cen. We searched among our retrograde 
stars for a possible identification of 10 Cen candidates (e.g. in 
the metallicity space) without any success. Therefore we favour 
small-number statistics to explain these large velocity disper- 
sions, although we cannot exclude their presence. 



7. Derivation of tlie radial scale lengths and scale 
heights 

Supposing that the thick disc and the thin disc are in equilib- 
rium, the velocity ellipsoids that were derived in the previous 
sections can be used with the Jeans equation to infer an estimate 
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Fig. 13. Eccentricity distributions for the thin disc (in red), thick 
disc (in dotted green) and halo (in dotted-dashed blue). The can- 
didates were selected based on their kinematics, according to 
the probabilistic approach of lRuchti"et al. (2010). This selection 
therefore introduces biases in the eccentricity distributions be- 
cause the tails of the velocity distributions cannot be selected. 



of their radial scale lengths {Iir) and scale heights (hz). Below we 
considered as thin disc the stars below 800 pc from the galactic 
plane, and as thick disc the stars between 1 and 3 kpc to avoid 
a strong contamination from the other components. In addition, 
the dispersions of the velocity ellipsoids were corrected from 
the observational error s by taking them out quadratically, as in 
I Jones & Walkeii (Il988l) . In cylindrical coordinates the radial and 
azimuthal components of the Jeans equation are 



pKz 



dlnR 



da-l 



dZ 



dZ 



1 dRpo-y^^ 
^ R dR 



(12) 
(13) 



where p is the density of the considered galactic component, 
Vc - 220 km s"' is the circular velocity at the solar radius, is 
the mean rotational velocity of the stars having the crvg, cry^,crv^ 

velocity dispersions, cr^^^ = VrVz-Vr Vz, and Kz is the vertical 



galactic acceleration. 



7.1. Radial scale lengths 
We consider that p{R) 



exp(-R/hii), and that cr^^ has 



the same radial dependence as p (as in ICaroUo et al.l l2010l) . 
Therefore, cr^^ oc exp( -/?//;«). In addition, one can assume that 
dcTy^ ^ /5Z w 0, which is true if the galactic potential is an infinite 
constant surface density sheet (iGilmore et al.ll989l) . In this case. 
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Table 7. Mean metallicity and kinematic values (not coiTected for observational errors) of the different galactic components (see 
Sect. l6.2l l in cylindrical coordinates. 



Galactic component 


N 


Vr 




Vz 


wim 












C [M/H] 






(kms-') 


(km s-') 


(kms-') 


(dex) 


(km s-') 


(km s 


-') 


(km s 


-') 


(dex) 


Thin discj.,,,5 


154 


-8 ±2 


-211 ± 1 


-5 ± 1 


-0.22 ± 0.02 


38 + 2 


26 + 


1 


20 + 


1 


0.28 + 0.02 


Thick discj.;„f 


193 


8±2 


-166 + 2 


-3 ±2 


-0.41+0.02 


56 + 4 


42 + 


3 


55 + 


2 


0.34 + 0.02 


Inner halo;;,,,^ 


105 


38± 13 


-60 ± 10 


-23 ± 10 


-0.70 + 0.03 


194 + 14 


121 + 


13 


122 + 


13 


0.59 + 0.05 


Thin discz 


163 


-1 ± 1 


-204 ± 1 


-2± 1 


-0.27 ± 0.02 


43 + 2 


33 + 


1 


25 + 


1 


0.31+0.02 


Thick discz 


187 


-2 + 3 


-167 ±3 


-12±2 


-0.45 ± 0.02 


66 + 5 


57 + 


4 


53 + 


3 


0.36 + 0.02 


Inner haloz 


45 


11 ±22 


-47 ± 21 


-48 ± 17 


-0.92 + 0.06 


215 + 31 


166 + 


28 


146 + 


22 


0.76 + 0.06 



Notes. Thin discz, thick discz and inner haloz are defined as the stars lying below Z < 800 pc, between 1 < Z < 4 kpc and above Z > 5 kpc, 
respectively. 

Table 8. Same as Table|2]but in cartesian coordinates. 



Galactic component 


U 




V 




W 




o-u 




CTv 




aw 






(km s" 


') 


(km s" 


') 


(km s" 


-') 


(km s 


-') 


(km s 


-') 


(km s 


-') 


Thin disct,„(. 


-18 + 


2 


-14 + 


1 


-5 + 


1 


38 + 


2 


25 + 


2 


20 + 


1 


Thick disct,„j 


-40 + 


3 


-63 + 


2 


-3 + 


2 


58 + 


4 


40 + 


3 


55 + 


2 


Inner haloj.;,,^ 


-65 + 


13 


-175 d 


:9 


-23 + 


10 


208 + 


14 


97 + 


12 


122 + 


13 


Thin discz 


-20 + 


1 


-21 + 


1 


-2 + 


1 


43 + 


2 


32 + 


1 


25 + 


1 


Thick discz 


-34 + 


4 


-69 + 


3 


-14 + 


2 


70 + 


6 


55 + 


4 


51 + 


3 


Inner haloz 


-33 + 


29 


-191 + 


25 


-49 + 


38 


223 + 


39 


142 + 


28 


158 + 


34 



Notes. Thin discz, thick discz and inner haloz are defined as in Table]?] The mean velocities are given without taking into account the solar motions 
(f/o, V^, Wo). 





Fig. 15. Metallicity and velocity distributions for the galactic components selected according to their position above the galactic 
plane. The Gaussian fits were obtained taking into account the errors on each parameter with 5 ■ 10^ Monte-Carlo realisations. 
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Fig. 14. Eccentricity distributions for the thin disc (in red), thick 
disc (in dotted green) and halo (in dotted-dashed blue). The can- 
didates were selected based on their distance above the galactic 
plane. 



the axes of the velocity ellipsoid are aligned with the cylindrical 
coordinates, and Eq.[T2lbecomes: 



It 



- 1 + 



2R 
'hR 



= 0. 



(14) 



An alternative option is to consider that the galactic potential 
is dominated by a centrally concentrated mass distribution and 
that the lo cal velocity ellipsoid points towards the galactic centre 
dGilmore et al 1989: Siebert et al..,2008.) . In that case, the previ- 
ous term becomes 



dZ 



1 - 



Equation [12] can then be rewritten as follows: 



■2-H 



2R 



Vr 



cr- 



0. 



(15) 



(16) 



Each of the terms of Eq. [T4]or Eq. [T6]were measured in our 
data, leaving as the only free variable the radial scale length 
liR of the discs. With the values derived for the thick disc of 
(cry„;cry^;\^) = (66 + 5;57 + 4; -167 + 3) kms"', we find 
Hr - 3.6 + 0.8 kpc using Eq. [141 and Hr ~ 3.4 + 0.7 kpc us- 
ing Eq. [161 These two values reasonably agree between them, 
and are found in the upper e nd of the values cit ed in the liter- 
ature (ranging fro m 2.2 kpc dCarollo et al.ll201C ) up to 3.6 kpc 
('Juric et al] |2008h . or even 4.5 kpc in the case of Chiba & BeersI 
(2001)). _ 

Using Eq. [161 and (o"v„; c^y,; V^) = (43 + 2; 33 ± 1; -204 + 
1) km s"', we find that the thin disc has a similar radial extent 
within our uncertainties as the thick disc, with /ir - 2.9 ±0.2 kpc. 
A smalle r thin disc has bee n suggested by other recent observa- 
tions (see lJuric et al.ll2008h . but once more, the value we derived 



is at the upper end of the previously reported values in the litera- 
ture. Still, an extended thin disc like this is plausible because our 
data mainly probe the old thin disc, which is likely to be more 
extended than its younger counterpart. 

7.2. Scale heights 

We assume that the last term of Eq.[T3lis negligible, because we 
are far from the galactic centre, and that p(Z) oc exp{-Z/hz). 
Equation[T3lhence becomes 



diner}, 
dZ 



1 Kz 
-H ^ = 0. 



(17) 



We us ed Kz = 2nG X 71 M^pc'^ derived bv lKuiiken & Gilmor j 
(|l99lh at |Z| - 1.1 kpc, but we note that this value of Kz 
might differ at the distances where our targets are observed. We 
also used for the thick disc the value derived from our data of 
dcTy^jdZ =15+7 km s"' kpc"' and cry^ =53 + 3 kms"'. 
Hence, for the thick disc, we find hz ~ 694 + 45 pc. 

We found for the thin disc that dcry^/dZ = 19 + 
10 km s"' kpc"' and cry^ =25 + 1 km s"', resulting in hz = 
216± 13pc. 

The derived values for both components agree well with, for 
example, lJuric et aL who suggested a thin disc with hz - 

300 pc, and a thick disc with hz = 900 pc. 



7.3. Metallicity dependence of the scale lengths of the thick 
disc 

To investigate furthermore the metallicity gradients found for the 
rotational velocity of the thick disc, we computed the radial scale 
lengths and scale heights for different metallicity bins using 
Eq. [island Eq. [T71 The results are shown in Table [9l where the 
metallicity bins were selected to include at least 30 stars each. 
Though we found that both /ir and hz increased with decreas- 
ing metallicity (except for the most metal-poor bin), this trend 
is not strong enough to stand-out significantly from the errors. 
We conclude that within the errors, the same scale lengths and 
scale heights are found, which is the signature of only one pop- 
ulation. Indeed, if a significant amount of relics of a destroyed 
m assive satellite s hould exist in our line-of-sight, as suggested 
by iGilmore et al.l ( |2002|) . one would expect them to have a dif- 
ferent spatial distribution compared to the canonical thick disc, 
which we do not observe. Unless, of course, the satellite debris 
provides the dominant stellar population in the thick disc. 

This result can also be discussed in the frame of a thick disc 
formed according to a radial migration scenario. In that case, 
the older stars that are at the solar radius have come from the 
inner parts of the Galaxy, and are expected to have a higher ver- 
tical velocity dispersion and a different metallicity, and hence, 
should exhibit scale heights dependent o n met allicity. In partic- 
ular, the model of Schonrich & Binne^ (120091) predicts a lower 
scale height for the metal-poor thick disc compared to its metal- 
rich counter part. This trend is not observed in our data (if it 
exists, it should be fairly small), which challenges the migration 
scenario as being the most important process of creation of the 
galactic thick disc. 



8. Conclusions 

A significant sample of roughly 700 low-resolution spectra 
(FLAMES/GIRAFFE LR8 setup, covering the IR Ca ii triplet) of 
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Table 9. Kinematic parameters, radial scale lengths and scale heights for different metallicity bins of the thick disc targets. 





N 


Vr 




V0 






Vz 




hR 


hz 


(dex) 




(1cm s-') 


(km s-') 


(km s 


-') 


(km s-') 


(km s-') 


(kms"') 


(kpc) 


(pc) 


-1.14 


36 


-5±9 


58+ 11 


-137± 


11 


61+ 7 


-7± 8 


59±7 


1.9± 0.7 


934± 166 


-0.67 


26 


-3± 17 


85± 17 


-161± 


11 


54± 11 


-4± 12 


54± 8 


4.0± 1.3 


804± 181 


-0.40 


56 


5± 8 


81± 8 


-168± 


6 


52±5 


-17± 6 


45±4 


3.8± 0.9 


610± 90 


-0.11 


37 


6±9 


64± 8 


-171± 


7 


50±6 


-18± 7 


45±5 


3.1±0.9 


620± 97 



Notes. All velocity dispersions were corrected for the observational errors. The scale lengths and scale heights were computed using 
Eq.[T6] and (TT] and we assumed Kz = 2nG x 71 Mgpc"^. 



stars faint enough to probe long distances and bright enough to 
get high signal-to-noise ratios were collected towards the galac- 
tic coordinates I ~ 277°, b ~ 47°. The stellar atmospheric pa- 
rameters (Tes, log g, overall metallicity [M/H]) were extracted 
from the spectra with our autom ated code, which is fully de- 
scribed in 



lOiha et al 



Kordopatis et alJ (1201 li) . Given the proper motions of 
(Il996h and our derived radial velocities, we were able 



to derive the distances and the positions for 479 stars of our sam- 
ple and the 3D motions and the orbital eccentricities for 452 of 
them. 

We found a thick disc with a mean rotational velocity ~ 
- 167 km s"' , a value slightly higher than the commonly adopted 
lag of less than 50 km s"' . We emphasise that our sample of stars 
probes distances above 1 kpc from the plane, so this difference 
may imply a correlation between lag velocity and vertical veloc- 
ity rather than a simple inconsistency with local data. The mean 
measured metallicity of the thick disc is -0.45 dex. Its metal- 
poor tails extends to - 1 .8 dex, whereas its metal-rich tail goes up 
to solar and super-solar values. Depending on how the thick disc 
stars are selected (Z positions or velocities), these values may 
vary a little, but generally agree well between them. The vertical 
velocity gradient, dV^/dZ =19 + 8 km s kpc and metallic- 
ity gradient <9[M///]/5Z = -0.14 + 0.05 dexkpc that are mea- 
sured for the regions where the thick disc is the dominant popu- 
lation (1 < Z < 4 kpc) seem to agree well with a smooth transi- 
tion between the ga lactic components, and are compatible with 



the valu es found bvlChiba & BeersI (l2()()ll):lA.llende Prieto et al 
(2006); Girard et al.1 (l2006h : iLee et al.l (l201 ih and iRuchti et~ 
( 1201 1). whose values range from 16 + 4 to 30 + 3 km s"' kpc"' 
and from -0.07 + 0.01 to -0.09 + 0.05 dex kpc"'. A corre- 
lation between the metallicity and the rotational velocity of 
dVild\MI H] = 43 + 11 km s '' d ex"' is also found, in agree- 
ment with lSpagna et al.1 (l20Toh and lLee eT al. (2011). 

The radial scale length and scale height of the thick disc are 
estimated to be /;« ~ 3.4 + 0.7 kpc a nd hz ~ 694 + 45 pc, which 
agrees well with the SDSS results of lJurii et an(l2008h . No clear 
metallicity dependences are detected for Hr or for hz, pointing 
towards a thick disc that is mainly composed by only one popu- 
lation. 

Finally, we found a broad peak of the eccentricity distribu- 
tion for t he thick disc arou nd e ~ 0.3 which, accordin g to the 
works of iSales et al.1 (l2009h and lDi Matteo etal.1 (l201lh seem to 
rule out a pure accretion sce nario. These res ults a gree with those 
rece ntly obtained by Wils on et al.l (1201 Oh and iDierickx et all 
(|2m0), who measured the eccentricity distributions for a sam- 
ple of RAVE and SEGUE stars. 

However, several questions still remain open. For instance, 
we found difficulties in fitting the transition between the thin and 
the thick disc simply by adjusting parameters in the Besangon 
model. In addition, the plateau for the high-metallicity and low- 
altitude stars (i.e. the thin disc) seems to suggest a local density 



of the thick disc around 18%, higher than that assumed by that 
model. Finally, the existence of intrinsic vertical gradients in the 
thick disc cannot be ruled out, because we did not obtain a suf- 
ficiently well defined (and statistically large enough) sample of 
thick disc members. Additional targets to increase the statistics 
and higher resolution spectra on well selected spectral domains 
to separate thin disc stars from the thick disc ones with chemi- 
cal content criteria are hence strongly recommended for future 
surveys. 
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Table 1. Kinematics of the selected targets belonging to the final catalogue 



ID 




Alt 






U 


Ac 


V 


A„ 


W 


Aw 


Vr 




V4, 




6 


Ae 




C/cen) 


C/cen) 


(km/s) 


(km/s) 


(km/s) 


(km/s) 


(km/s) 


(km/s) 


(km/s) 


(km/s) 


(km/s) 


(km/s) 


(km/s) 


(km/s) 






1 


0.142 


-0.361 


110.2 


3.8 


53 


50 


-125 


41 


27 


39 


-90 


50 


-66 


55 


0.54 


0.16 


2 


0.382 


0.107 


141.3 


5.3 


132 


68 


-47 


48 


127 


43 


-212 


69 


-73 


56 


0.59 


0.14 


3 


0.361 


-0.263 


9.0 


4.5 


8 


3 


-8 


4 


3 


4 


-25 


3 


-216 


4 


0.06 


0.01 


4 


0.354 


-0.739 


24.7 


6.4 


114 


51 


-134 


42 


-106 


39 


-149 


49 


-28 


48 


0.44 


0.08 



Notes. yU; and /ij are taken from lOiha et al.l ( 119961) . U, V and W are the cartesian velocity coordinates, with respect to the LSR, hence, without 
taking into account the peculiar solar velocities. Vr and are the velocity components in cylindrical coordinates, centred at the galactic centre. 

Table 2. Individual atmospheric parameters (as derived by the procedure of Paper I), relative V-Magnitudes, (B-V) colour, absolute 
V-magnitude and estimated S/N of the observed targets of the final catalogue 



ID 




log 8 


Alog g 


[M/H] A[M/H] 


niv 


(B-V) 


Mv 


Am,, 


S/N 




(K) 


(K) (cm s-2) 


(cm s"^) 


(dex) (dex) 










(pixer') 


1 


5962 


212 4.02 


0.30 


-0.76 0.15 


17.37 


0.61 


3.98 


0.86 


23 


2 


5001 


74 3.51 


0.12 


-0.25 0.10 


17.24 


0.69 


3.01 


0.41 


29 


3 


4272 


64 4.86 


0.12 


-0.22 0.09 


16.41 


1.42 


8.30 


0.13 


60 


4 


4997 


88 3.73 


0.17 


-1.00 0.16 


16.81 


0.73 


3.14 


0.29 


35 


Notes. 


m,. and fB-V) are taken fromlOiha et 


all (119961). 












Table 4. Positions of the selected targets belonging to the final catalogue 








ID 


/ 


b D 


Ad 


X Ax Y 


Ay 


Z 


Az 


R 


Ar 




n 


n (pc) 


(pc) 


(pc) (pc) (pc) 


(pc) 


(pc) 


(pc) 


(pc) 


(pc) 


1 


280.201 


47.452 4573 


2061 


552 245 -3069 


1365 


3397 


1511 


8167 


326 


2 


280.338 


47.363 6725 


1333 


816 159 -4477 


876 


4942 


967 


8504 


329 


3 


280.308 


47.406 401 


24 


48 3 -267 


16 


295 


18 


7955 


2 


4 


280.223 


47.409 5202 


705 


628 84 -3486 


467 


3853 


516 


8167 


124 



Notes. / and b are the galactic coordinates taken from lOiha et al.l ( 119961) . D is the line-of-sight distance. X, Fand Z are heliocentric distances in a 
cartesian reference frame. R is the galactocentric planar radial coordinate. 



